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^.^ ' Abstract 

i^H ' A phenomenological formula based on liquid drop model has been 

I , proposed for ground state binding energies of nuclei. The effect due to 

fj ■ bunching of single particle levels has been incorporated through a term 

^ ' resembling the one-body Hamiltonian. The effect of n-p interaction has 

P5 , been included through a function of valence nucleons. A total of 50 pa- 

' rameters has been used in the present calculation. The r.m.s. deviation 

. I for the binding energy values for 2140 nuclei comes out to be 0.376 MeV, 

►^ , and that for 1091 alpha decay energies is 0.284 MeV. The correspondence 

^^ ' with the conventional liquid drop model is discussed. 

lO ■ 

T— I ■ 1 Introduction 

t^ ■ 

f~^ ■ The charged liquid drop model of the nucleus, in combination with various 

f^ I modifications, has very successfully been applied to describe nuclear ground 

^-H . state binding energy and many other ground state properties [Ij. Most of the 

^ I modern global mass formulas retain the original Bethe-Weizsacker (BW) form 

as basis and add certain effects such as shell correction, deformation, surface 
symmetry terms, etc. Fully microscopic approaches such as Skyrme Hartree 
j^ ■ Fock or Relativistic Mean Field methods also have been employed to estimate 

the mass values through out the periodic table. Local formulas, as the name 
suggests, predict the unknown mass of a nucleus by extrapolation from the 
known masses of neighbouring nuclei. 

Fully microscopic approaches are yet to reach the accuracy achieved by 
microscopic-macroscopic mass formula. Skyrme Hartree Fock mass values show 
an r.m.s. deviation slightly less than 0.6 MeV[5]. Local approaches, on the 
other hand, can be applied to unknown mass regions only with limited con- 
fidence. Thus, for example, the neutron rich regions, important for various 
astrophysical processes such as r- or s-process, are usually studied using the 
microscopic-macroscopic approaches . 

Microscopic-macroscopic methods usually use the liquid drop model (LDM) 



X 



as the macroscopic part of the prescription and introduce different corrections 
based on microscopic models. Myers and Swiatecki[3] included the shell effects 
by considering the bunching of the single particle levels. Strutinsky correction [J 
takes care of the shell effects in a fully microscopic way. Finite Range Droplet 
Model (FRDM)[5] is an approach which, besides including Strutinsky correc- 
tions, pairing and other effects, also modifies the basic liquid drop macroscopic 
part of the mass formula. Another very successful approach based on a micro- 
scopic picture by Duflo and Zuker[B| yields an r.m.s. deviation of 0.373 MeV 
for 2135 nuclei. 

Estimation of binding energy remains a very important task in nuclear 
physics. For example, one of the most important inputs in the study of as- 
trophysical reactions is the nuclear ground state binding energy. In nuclei far 
from the stability valley, experimental mass measurements are very scarce and 
likely to remain so in near future. Thus, theoretical prediction of mass remains 
our only guide. In the present work we present a phenomenological formula for 
calculation of binding energy values. 

2 Theory and results 

2.1 The binding energy formula 

The general approach to a successful nuclear mass formula has been to assume 
the liquid drop model and add shell corrections to it. The liquid drop model 
gives a smooth curve while the shell corrections represent the effect of the single 
particle states. In the present work, we have employed a approach which is 
close to the Myers and Swiatecki prescription [31 with the shell correction being 
simulated using a number of parameters. 

The original liquid drop model mass formula is known to give incorrect re- 
sults because of the bunching of the single particle levels. To incorporate this 
effect, we have added a term which resembles the form of the one-body Hamil- 
tonian. However, the single particle energies are not constants throughout the 
periodic table but vary with the neutron and proton numbers of the concerned 
nucleus. To simulate this effect, we have used a term of the form 
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where i = 1,2 refer to the neutron or proton, respectively. The number of 
neutrons and protons in the nucleus are given by A/'^(= N) and Af'^{— Z), 
respectively. We assume a shell model like filling of the orbits to obtain the n^ 
values for different nuclei. Thus we have 



K+i~K 


for 


AA^ > K+i 


A/"' - Af^ 


for 


K<-^'<K+i 





for 


M' <MX 



Table 1: Values of various parameters for the binding energy formula obtained 
in the present calculation. The parameters kv,ks,kw are diniensionless while 
To is given in fm. The rest of the parameters are in MeV. 
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J^l 


ei 


AA^ 


4 


Uv 


11.8906 


1 


8 


0.1284 


8 


0.2027 


asurf 


7.3654 


2 


14 


-0.0227 


14 


0.0284 


fxy 


1.6627 


3 


20 


0.0651 


20 


0.0889 


ks 


-2.6383 


4 


24 


0.0306 


24 


0.0462 


w 


35.6813 


5 


28 


0.0394 


28 


0.0307 


kw 


8.3735 


6 


32 


0.0221 


40 


0.0191 


O-s 


3.1619x10-8 


7 


40 


0.0102 


44 


0.0256 


ro 


1.0625 


8 


44 


0.0248 


50 


0.0120 


< 


-6.8171x10-3 


9 


50 


0.0089 


64 


0.0106 


a-n 


-1.0645 


10 


60 


0.0138 


74 


0.0079 


ap 


-0.9508 


11 


64 


0.0093 


80 


0.0050 


alp 


3.3594x10-2 


12 


76 


0.0118 


82 


-0.0001 


a3 


-3.2543x10-"* 


13 


82 


0.0031 


88 


0.0013 


as 


4.7686 


14 


88 


0.0079 


96 


-0.0022 


a's 


-1.4162 


15 
16 

17 
18 


96 
122 
126 
132 


0.0074 
0.0041 
0.0043 
0.0078 







ei= 



e^= 0.0197 MeV 
0.0534 MeV and el= 



for Z < 34; iV > 46 
: 0.0312 for Z< 22;A^ 



> 29 



The different J\f^ values are estimated in the following way. While fitting the 
different parameters in the formula for binding energy, we have checked the 
errors at different neutron and proton numbers i.e. for different nuclei. When- 
ever, the errors tend to be large at a particular N or Z value, we assume that 
the single particle structure has changed at that particular number and new 
energy level has been occupied. Thus we have assumed Cq, values at different 
neutron or proton numbers (A/"^) as given in the last five columns of table 1. 
This particular form remains an ad hoc assumption in our work; yet we find 
that it can successfully explain the binding energy values. 

Some of the J\f^ values are easy to understand. The conventional magic 
numbers and semi-magic numbers do occur in the scheme naturally. On the 
other hand, some other values may be associated with subshell closures. For 
example, A/'q, = 14, 32, 64 and 76 (some in the case of neutrons only) agree 
with the 1^5/2, lp3/2, 2^5/2 and I/111/2 subshell closures. Similarly, A^ = 88 
may be associated with a deformed shell closure. Weaker deformed shell closure 
observed at A^ = 98 and Z = 100 are also very close to the value Na = 96 that 
we have obtained. The readers are referred to Refs.[7l [8j |9] for more details 



on shell closures and deformed magic numbers. We should also mention that 
Schwierz et al. [10] in a study on Woods Saxon potential have found indications 
of smaller gaps a.t N, Z = 14 and 32 in the light mass region. Sorlin et a/.|ll) 
have pointed out the evidence for subshell closures at 14, 32 and 64. 

We have also introduced modifications in a few e values and the A/"^ values 
at which they occur for a number of nuclei. This is done by again looking at 
the deviations and checking whether there is any systematic pattern. In nuclei 
with Z < 22 and iV > 29, we find that a better result is obtained if we change 
the €2 and eg values, i.e. for proton numbers corresponding to 14 and 20, 
respectively. Similarly, in nuclei with Z < 35 and N > 46, we find that a better 
result is obtained if we change the value of the parameter e| (corresponding to 
neutron number 44) and assume that the next change in single particle level does 
not correspond to A^ = 50. Interestingly, this may correspond to the possible 
disappearance of the magic gap at TV = 50 in nuclei with Z < 35. A similar 
situation appears in nuclei with Z < 42 where the bunching at A^ = 60 seems 
to disappear. The fact that all the above changes occur in neutron rich nuclei is 
a signature of rearrangement of single particle levels expected in such regions. 

The Binding energy formula used in the present work is given by 

B.E.{N, Z) = BldM + Bbunc + Bnp + Bel (2) 

BlDM — ai;(l — kyl )A — asur/(l — ksI'')A ' — Bsym — Bcoul 

^pair \'^) 



Bsym = W\I\{1 ~ kw\I\) + as{N- Zf (4) 

_3 Z(Z-iy ,.^^5/2 ..^ 

Bparr = ^.(^^^ + -^ + a,—^ (6) 

B„p = a^Np + ap7V„ + a\^{Np + N,,f + a\^{Np + N,,f (7) 

Here, (5„(p) is or 1, according as N{Z) is even or odd. We write / = 
(N — Z)/A. The expression for Bfj^nc is already given in eqn (1). In B„p, Np 
and Nn refers to the number of valence protons and neutrons. The electronic 
binding energy is estimated by the empirical relation 

Bei(MeV) = 1.44381 x lO^^Z^'^^ + 1.55468 x lO'^^Z^'^^ (8) 

The liquid drop part B^dm contains a number of correction terms along 
with the original BW prescription. As B},unc contain a large number of e^'s, we 
have a large number of parameters, viz. fifty in number. If we also consider that 
the nucleon numbers for the bunching of the levels as parameters, the number of 
increases substantially. One obvious problem in our approach is that the nucleon 
numbers used in the bunching part of the formula do not come naturally but have 
to be introduced from experimental mass values. Only microscopic calculations 
can possibly supplement experimental inputs in this regard. 



The bunching term aheady includes a significant part of the total energy. 
Thus we expect a major modification of the liquid drop model coefficients as 
well as the functional form of some of the smaller correction terms. We have 
included some different functional forms for the latter terms as described below. 

The volume and the surface terms as well the contribution of symmetry 
energy to them are retained in the LDM part of the formula. However, the 
parameters are expected to be substantially modified in the present approach. 
Another major modification that we find is in the symmetry energy term Bgym- 
The term proportional to |/| is the Wigner term that was introduced from the 
symmetry argument. We have introduced two additional terms in the symmetry 
as is clear from the expression. FRDM suggests that the charge asymmetry 
contribution should be large in heavier nuclei and uses a term proportional to 
{N — Z). We find that a term (iV — Z)^ fits the data better in our case. 

The Coulomb energy is represented as a sum of two terms. The first one 
is the usual Coulomb term for a charged liquid drop. For the second term, 
we have used a functional of Z"^ j A. The usual Coulomb contributions include 
an exchange term of the form Z^^'^ /A^/^, proton form factor correction of the 
form Z"^ I A. The FRDM also has contributions from volume and surface rear- 
rangement terms which depend on Z. Lunney et aZ. [T] have also used a term 
proportional to Z"^ j A? in the macroscopic part of the energy functional related 
to surface stiffness. We have included a term of the form iZP' j A^y'l'^ as a single 
term representing the Coulomb corrections. We should mention that a nearly 
comparable fit can be achieved using a form (Z^/A)^ also. 

Another addition is the Bn-p term which is a part of the microscopic en- 
ergy. It has been demonstrated [12j [13] that the number of valence neutrons and 
protons show a decidedly strong correlation with n-p interaction and thus with 
deformation. We have taken a form of the type shown in eqn (7). The number 
of valence particles or holes has been calculated from the nearest major shell, 
i.e. 8, 20, 28, 50 82, 126 and 184. We have also tried to fit the data using 172 
as the last major shell. We find that the fit is slightly worse in this case. How- 
ever, in keeping with the expectation that the shell closures are less pronounced 
in the SHE region, the binding energy data is not really sufficient to conclude 
whether the next major closure is at A^ = 184 or 172. On the other hand, as 
explained later, we have also calculated the Q- value in a-decay where the former 
value gives decidedly better results and we have used it in our calculation. The 
dependence of nuclear mass on valence nucleons has already been studied by 
Mendoja-Temis et a/. [14] in a different form. 

The pairing energy has been shifted so that the value for Bpair is zero in 
even-even nuclei. We have also assumed a A^^^^ dependence for the n-p pairing 
term. 

2.2 Parameters 

The justification of using such a large number of parameters can only be through 
a good and simple prediction. One advantage of using this formula is that one 
may use it directly without microscopic correction as the bunching due to shell 
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Figure 1: Deviations of binding energy for 2140 nuclei with N, Z > 8 from the 
predictions of the present formula. 



effects are included in the Bbunc term. We have fitted the experimental masses 
in two steps as explained below. 

In the first step we have included all the experimental mass values in nuclei 
with N,Z >8 from the mass table by Audi et aL[T5]. A total of 2140 nuclei have 
been used in the fitting procedure yielding an r.m.s. deviation of 0.378 MeV. 
More importantly, there are only 39 nuclei which show deviation larger than 1 
MeV and none more than 2 MeV. In nuclei with mass more than 100, only 4 
nuclei show a deviation larger than 1 MeV and none more than 1.5 MeV. We 
would like to note that if the modifications in the three e^ values, as explained 
in the previous subsection, are not introduced for neutron rich nuclei, the r.m.s. 
deviation becomes 0.400 MeV. In this case we also assume that for the nuclei 
with N < 22 and 34, the bunching at iV = 28 and 44, respectively, disappears. 

It should be interesting to compare the results for some measurements per- 
formed later than 2003 and not included in the mass table[r5] mentioned above. 
The masses of a number of nuclei have been measured either for the first time or 
with increased precision in the last few years. We have looked for measurements 
which are completely new or which substantially differ from previous values at 
least in some of the nuclei studied. All the measurements from such references 
have been included for completeness. Refs. [T1]~[3T] have been selected for in- 
clusion. These measurements have been divided in two groups. The first group 
contains those nuclei for which experimental values were available in Audi et 
a/.|15) but later measurements have modified them. There are 134 nuclei in 
this group. In the second group, there are 43 nuclei whose experimental bind- 
ing energy values have been measured later and were not tabulated in Audi et 
aZ.[T5]. For both the groups taken together, the r.m.s. error is 0.438 MeV. For 
the second group alone, the r.m.s. deviation is 0.517 MeV. 

In the second step, we have replaced the binding energy values of the nuclei 
belonging to the first group with the improved binding energy measurements 
obtained later wherever available. A fitting procedure for 2140 nuclei now show 
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Figure 2: Deviations of binding energy for 43 nuclei whose binding energy values 
have recently measured from the predictions of the present formula. 

Table 2: Deviation in some global mass formulas 







No 


of nuclei 


No. of 
parameters 


r.m.s. 
deviation 

(MeV) 


Myers- S wiatecki |35j 




1312 


50 


0.663 


Skyrme-HFB[2] 






2149 


24 


0.581 


Duflo-Zuker[B] 






2135 


28 


0.373 


FRDM[5] 






2135 


33 


0.676 


Present 






2140 


50 


0.376 



an r.m.s. deviation of 0.376 MeV. The quality of fit is as good as obtained in 
the first step. The new parameters are very close to the fitted values in the 
first case. The new fitted parameters also describe the values in the second 
group slightly better, the r.m.s. deviation being 0.494 MeV. Fig. 1 shows the 
deviation from the prediction of the (2) formula for the 2140 nuclei while Fig. 
2 is the corresponding figure for the 43 new measurements. The parameters are 
given in table 1. 

In table 2 we compare our results with some other global calculations. It 
is clear that our calculation uses more free parameters. Besides, a number 
of parameters in the other approaches are predetermined from considerations 
other than mass measurements. However, as already noted, the advantage of 
the present formula is that it is an analytic expression, and not dependent on 
input from any other calculation. 

To check the agreement of our results with other predictions, we have chosen 
two isotope chains, Z = 50 and Z = 110. The former has been chosen because 
experimental binding energy values for a long chain of isotopes, from A^ = 50 to 
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Figure 3: Differences of the present calculation from experimental values and 
other prescriptions for Z = 50 and Z = 110 nuclei. See text for details. 



TV = 84, are available. The other has been chosen to compare the predictions 
in a region where no experimental data is available. The percentage deviation 
values for binding energy are plotted in Fig. 3. The legends AW correspond to 
the experimental or estimated values from the mass table[T5]. The legends HFB, 
FRDM and DZ refer to the calculations according to the Skyrme Hartree Fock 
Bogoliubov calculation [2^, FRDM calculation j5j and Duflo Zuker ten parameter 
formula [6]. 

The figure indicates that the present formula agree well for the Z=50 iso- 
tope chain. However, the right panel show that there is a systematic deviation 
between the present result and the predictions of HFB, FRDM or Duflo-Zuker 
calculations for Z = 110 isotopes. The present model systematically predicts 
less binding energy for all the nuclei in the chain. More importantly, beyond 
N = 172, the deviations with all the three theoretical calculations start to rise. 
When experimental mass values for these nuclei become available, it will be 
clear whether the present prescription needs to be modified. It needs to be re- 
membered that the largest values that have been used in the fitting procedure 
are TV = 159 and Z = 108, as available in Audi et al.fTE\. 

As a further check on the present prescription, we have calculated the alpha 
decay energies in 1091 decays obtained from Audi et al.^5\ augmented by the 
latest results available in the National Nuclear Data Center website |32]. The 
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Figure 4: Error in alpha decay energy for 1091 alpha decays 



r.m.s. deviation for all these decay energies is 0.284 MeV, a value consistent 
with the accuracy in the prediction of binding energy in the present work. The 
deviations in alpha decay energies are plotted in Fig. 4. 

It is possible that beyond the highest nuclcon numbers with known exper- 
imental binding energy values, i.e. N — 159 and Z = 108, there are other 
nucleon numbers where we may have to introduce additional bunching. We 
have looked at the alpha decay energies for nuclei with larger N or Z values and 
have tried to fit with different e parameters at different N or Z numbers. A least 
square fitting did not yield a substantial improvement in the alpha decay ener- 
gies. However, as the number of observed decays in this region is not large, it 
is not possible to conclude positively about possible change in bunching beyond 
the above N or Z values. We have already noted that alpha decay Q-values 
indicate that A^ = 184 is a better closed core than 172. 

2.3 Correspondence with the hquid drop model 

We see that the usual liquid drop model parameters are modified to a large 
extent in the present calculation. It is difficult to compare the standard values of 
parameters with the specific form in (1) . Parts of all the conventional liquid drop 
model terms, such as the volume, surface, symmetry and Coulomb energies, are 
included in the bunching part of the formula. Thus the values of the different 
parameters in the present work differ substantially from the values obtained 
from their counterparts in conventional calculations. Even a parameter like rg, 
which seems to take a value in conformity with other results, is actually very 
much dependent on the form of the functional chosen for the second term in 

Bcoul- 

It is possible to visualize our results in the following way. In N — Z nuclei, 
there is no contribution from symmetry energy. We would like to interpret 
the difference in contribution from the proton and the neutrons to B^unc in 
such nuclei as essentially a part of the Coulomb energy in the conventional 



LDM. Thus we have written the Coulomb energy in TV = Z nuclei, neglecting 
the difference between the neutron and proton energies in Bnp, as Ecoui = 
Bcoui + B^^^^ — -Bfetmc- This quantity we have calculated for A^ = Z nuclei up 
to 108 from the proposed binding energy formula using the parameters deduced 
above and have plotted it with a simple conventional form consisting of a direct 
term and a diffuseness correction, i.e. 

-^^ + /- (9) 

Royer and Rousseau [33] pointed out that the Coulomb exchange term is unnec- 
essary if the diffuseness correction is included. Fitting Ecoui with the above 
expression, we obtain ro=1.2114 fm and /p = 1.0395 MeV which are within the 
range of usual values [33]. The Coulomb energy Ecoui is better described using 
a functional of the form 

3 e^Z^ ^ Z^ Z^ 

In even-even N ^ Z nuclei, it is possible to extract the volume and surface 
term coefficients from the mass formula. The sum of volume and surface terms is 
obtained by adding Ecoui to the total B.E. from eqn (2) taking the microscopic 
corrections into account and fitted with a functional of the form 

ay A - a^A2/3 (11) 

The values of the coefficients ay and as have been obtained in this way. We 
have considered only nuclei for which the microscopic correction terms from [5] 
are available. 

With the equivalent volume, surface and Coulomb energies given by the 
prescriptions (10) and (11), it is possible to estimate the symmetry energy con- 
tained in the values obtained from eqn (2). We have calculated the symmetry 
energy for all the even even nuclei among the 2140 which we have fitted origi- 
nally taking into account the microscopic corrections from Moller et al. [5] . There 
are various possible expressions for the symmetry energy in LDM and we have 
taken a functional of the form 

- aykyl^A + asksI^A^'^ - W\I\ + Cs{N - Z) (12) 

the last term being used in the FRDM. The different parameters have been 
estimated by least square fitting. 

The interpretation of pairing energy terms for odd-Z or odd-N terms are 
straightforward. The coefficient of the 1/ yjN{Z) term is 4.77 MeV, a value 
very close to the accepted value of 4.8 MeV. The n-p pairing term is usually 
taken in a different form from the present work and cannot be directly compared. 
Its effect is of course small. 

Following the procedure described above, we get the following expression for 
binding energy in terms of the deduced parameters: 

B(MeV) = 15.2893(1 - 1.8735/2)A - 16.4662(1 - 2.4595/2)^2/3 

10 



-foM4W + 0-3639^ -0.04211, 
-28.0438|/| - 0.4857(iV - Z) ~ Bparr - E.heii (13) 

where Esheii values are the microscopic corrections from Moller et al.^. It is 
seen that the coefficients in the above formula are within the accepted ranges 
for their counterparts in the liquid drop model. We note that the parameters in 
the above expression have been obtained fitting, not the experimental values, 
but the values predicted by eqn (2). Using the above formula, we find that 
the r.m.s. deviation from experimental values is obtained as 0.969 MeV for the 
nuclei in Audi et al.^^ Thus, we may conclude that the present prescription 
is equivalent to a conventional LDM formula which takes shell corrections into 
account. 



3 Summary 

A new phenomenological formula has been proposed for ground state bind- 
ing energies of nuclei. Two additional components have been included in the 
proposed formula in addition to the conventional liquid drop model. The con- 
tribution of the bunching of single particle levels has been included through a 
term resembling the one-body Hamiltonian. No additional shell correction from 
any external calculation has been considered. The effect of n-p interaction and 
deformation has been included through a function of the number of valence nu- 
cleons. A total of 50 parameters has been used in the present calculation. The 
r.m.s. deviation in the ground state binding energy for 2140 nuclei comes out to 
be 0.376 MeV. As an additional check on the formula, decay energies for 1091 
alpha decaying nuclei has been calculated which yields an r.m.s. deviation of 
0.284 MeV. As the bunching term incorporates a large part of the conventional 
liquid drop model terms, the parameters in the present work differ substantially 
form their standard values. The present formula may be reduced to the form of 
a conventional liquid drop model with parameters close to standard values. 
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